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ABSTRACT 

Above  ground  structures  which  have  been  officially  designated  as 
fallout  shelters  pose  a  number  of  ventilation  problems  which  require 
attention  to  insure  that  the  thermal  environment  of  the  protected  area 
will  be  habitable.  The  various  materials  and  configurations  of  the 
structures  and  the  effect  of  solar  radiation  requires  that  the  heat 
transfer  through  the  walls  and  other  surfaces  be  considered  separately. 
To  provide,  heat,  transfer  data  for  those  structures  which  are  of  thick 
wall  construction,  a  widely  accepted  analytical  solution  was  programmed 
for  the  1820  computer. 

i 

'  A  modified  psychrometric  chart  was  developed  so  the  sensible  heat 
factor  technique  can  be  used  to  determine  ventilation  requirements  for 
above  ground  shelters  subjected  to  unusual  climactic  conditions. 

Sample  calculations  for  a  500-man  shelter  located  in  St.  Louis, 
Missouri  show  that  the  maximum  heat  gain  through  the  thick  walls  is 
only  1.79  per  cent  of  the  human  load  and  the  heat  loss  through  the  floor 
is  3.33  per  cent  of  the  human  load.  The  people  in  this  case  contribute 
almost  the  entire  net  heat  load.  Calculations  also  indicate  that  the 
required  ventilation  rate  is  almost  twice  that  required  for  the  same 
shelter  in  Sacramento,  California.  This  difference  underscores  the  fact 
that  frequently  quoted  rates  from  3  to  18  cfm  per  person  can  be  mislead¬ 
ing  unless  properly  qualified. 

Continued  work  on  this  task  is  directed  toward  the  accumulation  of 
more  data  on  heat  transfer  through  walls  of  heavy  construction,  and  heat 
loss  through  the  floor.  Further  modifications  to  the  psychrometric 
chart  may  be  needed  and  the  inside  design  conditions  are  to  be  investi¬ 
gated  with  respect  to  comfort  vs,  economy. 
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INTRODUCTION 


Many  buildings  now  in  existence  have  been  designated  as  fallout 
shelters  by  Civil  Defense  authorities,  but  only  when  they  are  properly 
ventilated  can  they  realistically  be  classified  as  protective  areas. 

Most  of  the  designated  buildings  are  of  conventional  heavy  construction 
with  masonry  walls  12  to  16  inches  thick;  however,  in  the  event  of  a 
nuclear  threat  it  would  be  possible  to  thicken  many  of  the  walls  with 
sand  bags,  rows  of  concrete  block  or  extra  layers  of  concrete.  There* 
fore,  the  first  phase  of  this  study  on  ventilation  for  above  ground 
shelters  is  directed  toward  thick  composite  wall  construction.  Recogni¬ 
zing  the  likelihood  of  non-uniformity  in  the  walls  of  such  shelters, 
the  solution  treats  the  walls,  roof  and  floor  as  separate  components  of 
any  thickness  or  material  giving  complete  flexibility  to  the  solution. 

The  effective  temperature  in  a  protective  shelter  above  ground  is 
influenced  by  the  following  factors: 

a.  Internal  heat  load 

1.  Humans  are  the  prime  source  of  both  sensible  and  latent 
heat  inside  the  shelter.  Items  such  as  lights,  electronic  equipment 
and  fan  motors  produce  sensible  heat  only. 

b.  Thermal  characteristics  of  the  soil 

1.  Above  ground  shelters  having  floors  in  direct  contact  with 
the  ground  lose  heat  through  the  floor  to  the  soil.  The  rate  of  heat 
transfer  depends  on  the  thermal  characteristics  of  the  Boil  and  the 
difference  between  the  floor  and  soil  temperatures.  During  the  first 
few  hours  of  occupancy  the  rate  of  transfer  is  quite  high  but  it 
rapidly  decreases  and  approaches  a  constant  value  in  10  to  20  days. 

c.  Thermal  characteristics  of  the  building  materials 

1.  The  thermal  characteristics  of  the  building  materials  do 
not  vary  much  from  one  shelter  to  another  because  in  all  cases  the  mass 
of  the  materials  is  necessarily  high  for  protection  against  nuclear 
radiation. 

d.  Sol-Air  temperature 

1.  The  sol-air  temperature  can  be  readily  calculated  and  does 
not  vary  a  great  deal  from  one  location  to  another  within  the  latitudes 
bounding  the  Continental  United  States  of  America. 


e.  Condition  of  the  outside  air 


1.  The  design  condition  of  the  outside  air  to  be  used  for 
ventilation  is  based  on  climatological  data  gathered  over  a  long  period 
of  time  and  plays  an  important  part  in  controlling  the  effective  temp¬ 
erature  inside  the  shelter. 

The  factors  mentioned  here  involve  numerous  variables  and  it  should 
be  emphasized  that  the  problem  is  complex  and  it  is  highly  unlikely  that 
it  will  ever  be  resolved  into  a  simple  set  of  rules  or  tables  applicable 
to  all  shelters. 

For  this  report,  only  walls  of  3-foot  thickness  have  been  considered, 
but  the  analysis  can  be  applied  to  any  thickness  providing  there  are  no 
more  than  3  layers  of  materials. 


SOLUTION  DESCRIPTION 

The  ventilation  requirement  for  an  above  ground  structure  is  based 
on  the  heat  gain  through  the  exposed  areas  plus  the  internal  heat  load 
minus  the  heat  discharged  through  the  floor.  If  the  ventilating  air 
has  the  effect  of  adding  heat  rather  than  removing  heat  it  should  then 
be  kept  to  a  minimum  and  some  form  of  refrigeration  employed. 

Heat  Gain  Through  Exposed  Walls  and  Roofs 

Periodic  heat  flow  through  composite  yalls  and  roofs  has  been 
thoroughly  explored  by  Mackey  and  Wright,  and  £he  results  of  their 
work  form  the  basis  for  the  cooling  load  tables^  adopted  by  the  American 
Society  of  Heating,  Refrigerating,  and  Air  Conditioning  Engineers. 

These  tables  are  confined  to  conventional  walls  and  roofs;  consequently, 
for  thick  wall  construction  it  was  necessary  to  compute  the  cooling  load 
by  utilizing  the  mathematical  solution  of  Mackey  and  Wright.  In  this 
solution  the  harmonic  decrement  factor  (A)  and  the  harmonic  lag 
e  (  0  ),  (given  in  Appendix  A),  are  arithmetically  cumbersome, 

th  .t  these  parameters  were  programmed  for  the  1620  computer.  Equations 
to  determine  the  temperature  of  the  inside  wall  surface  (  T  )  and  the 
instantaneous  rate  of  heat  transfer  (  Q/A  )  (given  in  Appendix  A) 
were  also  included  in  the  computer  program. 

Complete  details  of  the  1620  computer  program  are  shown  in  Appendix 
A  making  it  possible  for  anyone  with  access  to  this  computer  to  deter¬ 
mine  the  periodic  heat  flow  for  any  wall  without  recourse  to  the  origi¬ 
nal  solution.  For  convenience,  24  walls  consisting  of  eight  cases  of 
concrete- earth  configurations  with  three  types  of  earth  for  each  case 
have  been  analyzed  for  two  inside  air  temperatures  of  85  F  and  90  F. 

The  eight  cases  are  shown  in  Figures  IB,  2B,  3B  and  4B  of  Appendix  B. 
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Cases  V  and  VI  would  most  likely  represent  roof  constructions.  Cases 
VII  and  VIII,  similar  to  V  and  VI,  might  represent  walls  where  the 
earth  is  in  the  form  of  sand  bags  or  heid  m  place  by  piywood  which 
would  have  1 i tt*i p  thermal  effect  on  s  thick  wall.  The  characteristics 
of  the  three  types  of  earth  (labelled  A,  B,  and  C)  are  given  in  Table  I. 
A  reference  to  Case  IA  means  Case  I  configuration  with  type  A  earth. 

The  thermal  characteristics  for  concrete  are  also  given  in  Table  I. 


Table  I.  Properties  of  Earth  and  Concrete 


Material 

Type 

k 

P 

C 

kpc 

Moisture 
(%’)  of 
dry  wt. 

Earth 

A-silty  clay 

.17 

75.3 

.193 

2.47 

5 

B-clay  sand 

.66 

107 

.163 

11.5 

5 

C-clay  sand 

1.74 

107 

.163 

30.4 

20 

Concrete 

Sand  and  gravel 
aggregate 

1.0 

140 

.25 

1 

1 

1 

1 

Cinder  aggregate 

.30 

100 

.20 

-- 

-- 

The  computer  produced  a  ponderous  set  of  results  for  the  two 
temperatures,  but  it  was  possible  to  group  the  results  into  the  two 
tables  shown  in  Appendix  C.  Using  these  tables  a  design  engineer  can 
determine  for  each  exposure  the  maximum  inside  wall  surface  temperature, 
the  maximum  heat  transfer  through  the  wallB  (positive  or  negative)  in 
BTU/hr/sq  ft,  and  the  corresponding  time  of  day  at  which  these  values 
occur.  Of  course,  his  parameters  must  satisfy  the  conditions  imposed 
on  the  development  of  the  tables  —  otherwise  the  computer  program  is 
required.  For  conventional  walls  the  tables  in  the  ASHRAE  guide  may 
be  used. 

Heat  Loss  Through  The  Floor 

For  a  plane  surface  in  soil  large  enough  so  that  edge  effects  can 
be  neglected  and  held  at  a  constant  temperature  differential  AT 
above  the  initial  soil  temperature,  the  instantaneous  rate  of  heat 
transfer  is  given  by  the  formula3 
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. 567  ykpc  AT 
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where  q  is  the  heat  transfer  from  one  side  of  the  plans  surface  in 
BTU/hr/ft^ 

k  =  thermal  conductivity  of  the  soil,  BTU/hr  ft  F 

p  ■  density  of  the  soil,  lb/ft3 

c  -  specific  heat  of  the  soil,  BTU/lb  F 

AT  =  temperature  difference  between  the  plane  surface  and 
the  initial  temperature  of  the  soil 

and 

t  ■  time  in  hours 

The  kpc  values  of  different  soils  are  readily  available  permit- 
ting  quick  comparisons  of  their  potential  heat  absorption  rates. 

Formula  (a)  is  somewhat  awkward  to' use  but  Raber,  BoeBter  and  ’ 

Hutchinson  have  plotted  it  for  convenience,  as  shown  in  Figure  1.  The 
use  of  the  formula  assumes  a  constant  AT  which  would  not  be  the  case 
during  the  initial  period  of  occupancy  when  the  floor  temperature  would  i 

be  rising.  However,  the  accuracy  of  the  method  increases  with  the  time 
of  occupancy,  so  for  periods  of  10  days  or  more  the  results  are  quite 
reasonable.  For  shelters  under  continual  use,  the  soil  may  have 
absorbed  sufficient  heat  to  nullify  its  effectiveness  as  a  heat  sink, 
thereby  obviating  the  necessity  for  this  calculation. 

Internal  Heat  Load 

People,  lights,  stoves  and  electronic  equipment  will  all  add  to 
the  internal  heat  load;  however,  under  the  austere  conditions  specified 
for  fallout  shelters  the  people  will  be  the  prime  contributors.  The 
sensible  and  latent  heat  emitted  by  the  human  body  varies  with  clothing, 
the  degree  of  physical  activity  and  the  wet  and  dry  bulb  temperatures. 

It  is  clear  that  body  heat  loss  cannot  be  accurately  forecast,  but  for 
illustration  purposes  the  sensible  heat  loss  will  be  estimated  at 
100  BTU/hr  and  the  latent  heat  loss  at  300  BTU/hr  assuming  people  at 
rest.  If  the  occupants  were  working  these  values  might  be  doubled  or 
tripled,  so  estimates  of  the  human  heat  load  must  be  made  with  some 
factor  of  safety.  The  heat  produced  by  auxiliary  equipment  will  be 
almost  entirely  sensible. 
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Net  Heat  Load 
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The  net  heat  which  must  be  removed  from  the  shelter  will  be  the 
algebraic  summation  of  heat  transferred  through  the  roof,  walls,  floor 
and  the  internal  heat  load.  The  net  amount  must  be  divided  into 
sensible  and  latent  heat  so  that  the  condition  and  amount  of  the 
ventilating  air  can  be  properly  determined.  If  outside  air  cannot  meet 
the  required  condition  then  no  amount  of  air,  within  practical  limits, 
can  cool  the  shelter  without  refrigeration.  The  erroneous  assumption 
is  sometimes  made  that  the  quantity  of  air  is  the  sole  criterion  in 
ventilation  design. 

Sensible  Heat  Factor 

The  conditions  which  will  most  likely  prevail  in  shelters  with 
thick  walls  result  in  low  sensible  heat  factors  which  can  not  be  handled 
graphically  on  all  psychrometric  charts.  To  overcome  this  problem,  a 
modified  psychrometric  chart  (Figure  2)  was  developed  to  include  low 
sensible  heat  factors  and  also  the  isotherms  of  several  effective 
temperatures.  Through  the  use  of  this  chart  the  condition  of  ventilat¬ 
ing  air  which  will  meet  the  heat  removal  demands  of  the  shelter  can  be 
graphically  determined  and  its  rate  of  flow  subsequently  calculated. 
Also,  the  refrigeration  requirements  can  be  determined. 


ILLUSTRATIONS 

To  clarify  this  procedure,  illustrations  are  presented  which  are 
concerned  with  two  identical  shelters  located  in  different  cities. 

Illustration  No.  1 


Consider  an  above  ground,  500-man  shelter  which  is  normally  unoc¬ 
cupied  and  located  in  St.  Louis,  Missouri.  The  inside  dimensions  are 
70-  by  70-  by  10- feet.  The  roof  is  Case  VA  and  the  walls  are  Case  IA. 
The  floor  is  a  concrete  slab  on  type  A  earth.  If  the  shelter  is 
occupied  by  500  persons  for  two  weeks  and  the  inside  conditions  are  not 
to  exceed  90  dry  bulb  and  85  effective  temperature,  how  much  ventilation 
is  required;  what  must  be  the  condition  of  the  air  entering  the  shelter; 
and  how  many  tons  of  refrigeration  are  required? 

Solution 


1.  The  first  step  is  to  establish  the  time  of  day  at  which  the  maximum 
cooling  load  occurs.  This  requires  a  trial  and  error  method  but  a  con¬ 
clusion  can  soon  be  reached  assuming  that  adequate  heat  transfer  data 
is  available.  For  example,  the  heat  gain  through  a  large  flat  roof 
will  usually  be  much  larger  than  that  through  the  walls  and  when  the 
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internal  load  is  essentially  constant,  as  in  illustration  No.  1,  the 
roof  becomes  the  determinant  factor.  By  referring  to  Table  IC  (Appendix 
G)  it  can  he  seen  that  the  maximum  roof  load  of  0.751  BTU/hr  sq  ft, 
resulting  from  a  peak  sol-air  temperature  in  the  early  afternoon  of  the 
previous  day,  will  occur  at  12  noon.  The  thick  wall  of  heavy  construc¬ 
tion  is  very  effective  in  modulating  the  cyclic  temperature  changes  of 
the  inside  wall  surface,  consequently  the  heat  transfer  rates  at  10  AM 
and  2  PM  are  practically  the  same  as  at  12  noon.  Note  that  the  north 
and  south  exposures  have  a  negative  heat  transfer  which  is  due  to  the 
high  inside  temperature  and  the  degree  of  solar  radiation  on  those 
surfaces.  If  the  12  noon  values  are  applied  to  the  areas  in  the  given 
problem  the  net  heat  gain  through  walls  and  roof  amounts  to  3579  BTU/hr. 


2.  The  second  step  is  to  determine  the  heat  loss  through  the  floor 
which  can  be  done  through  the  use  of  Figure  1  or  by  using  formula  (a). 
The  value  for  kpc  is  given  in  Table  I  and  the  time  given  is  336  hours 
(two  weeks) .  In  determining  a  value  for  AT  the  soil  temperature  was 
assumed*  to  be  57  F  which  is  the  well  water  temperature  in  the  St.  Louis 
area. ^  The  floor  surface  temperature,  which  is  originally  at  soil 
temperature,  rises  rapidly  after  occupancy  and  then  levels  off  in  an 
asymptotic  manner  as  It  approaches  the  room  air  temperature.  For 
practical  and  conservative  purposes  it  can  be  considered  a  constant 
value  approximately  five  degrees  below  the  room  air  temperature.  For 
this  illustration  the  floor  temperature  would  therefore  be  85  F  and 
AT  -  (85-57)  giving  a  total  heat  loss  to  the  soil  of  6660  BTU/hr. 


3.  The  third  step  is  estimating  the  internal  load  which  in  the  given 
illustration  involves  humans  only.  Assuming  the  occupants  are  Bitting 
or  reclining  in  light  clothing  the  sensible  load  has  been  estimated  at 
50,000  BTU/hr  and  the  latent  load  150,000  BTU/hr. 


4.  The  net  sensible  heat  load  is  therefore  3,579  -  6,660  +  50,000  *> 
46,919  BUT /hr  and  the  net  latent  heat  is  150,000  BTU/hr. 


5.  The  sensible  heat  factor  =  _ sensible  heat _ 

sensible  heat  +  latent  heat 


- 46jJ19 -  D  0.238 

46,919  +  150,000 


6  .  Reference  is  now  made  to  the  modified  psychrometric  chart  (Figure  2) 
and  the  following  steps  are  taken: 

a.  Notice  the  location  of  reference  point  "A"  which  is  an  integral 
part  in  the  construction  of  the  original  psychrometric  chart. 


This  assumption  is  based  on  information  given  in  the  ASHRAE  guide  for 
heating  load  calculations.  For  an  actual  shelter  local  soil  temperatures 
should  be  obtained. 
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b.  Locate  the  inside  design  conditions  of  90  dry  bulb  and  85 
effective  temperature  and  call  it  "B." 

c.  Locate  the  uuLeiue  design  temperature  for  St.  Louis,  94  dry 
bulb  and  77  wet  bulb  and  call  it  "C." 

d.  Draw  a  horizontal  line  from  "C"  to  the  saturation  line;  call 
it  "CD."  When  the  air  at  condition  "C"  is  coaled  it  will  move  left 
along  this  line  to  "D"  and  then  down  the  saturation  line. 

e.  Draw  a  line  from  "A"  to  the  sensible  heat  factor  value  of 
.238  and  call  it  "AE."  This  establishes  the  slope  which  the  air  will 
follow  as  it  heats  up  Inside  the  shelter. 

f.  Draw  a  line  through  "B,"  parallel  to  "AE"  cutting  "CD"  at  "F" 
which  represents  the  required  entering  condition. 

7.  When  the  air  at  condition  "C"  (St.  Louis  design)  is  cooled  either 
by  nature  or  refrigeration  it  must  reach  the  condition  at  point  "F"  in 
order  to  handle  the  internal  load.  The  entering  air  temperature  at  "F" 
has  a  dry  bulb  temperature  of  82.2  F  and  a  wet  bulb  temperature  of 
73.8  F. 

8.  The  quantity  of  air  required  is  calculated  by  the  following  formula 


sensible  heat 
crm  a  - 

1.08  x  ("B"  dry  bulb  -  "F"  dry  bulb) 

For  this  illustration  the  cfm  equals  5560  total  or  10.1  cfm/person, 

If  the  entering  air  temperature  moves  to  the  left  of  "F"  the  same  quan¬ 
tity  of  air  will  cause  the  inside  condition  to  move  to  the  left  of  "B" 
producing  a  more  comfortable  condition  in  the  shelter.  If  the  entering 
air  temperature  moves  to  the  right  of  "F"  refrigeration  is  needed. 

9.  The  refrigeration  required  to  cool  the  air  from  the, outside  design 
temperature  to  the  condition  at  "F"  can  be  computed  through  the  use  of 
enthalpies.  The  modified  psychrometric  chart  .shows  the  enthalpies  for 
points  "C"  and  "F"  to  be  40.2  BTU/lb  and  37.3  BTU/lb  respectively. 

The  refrigeration  required  is  therefore 

5560  x  .075  x  (40.2  -  37.3)  «  6.06  tons 

200 


Illustration  No.  2 


The  second  illustration  assumes  that  all  criteria  for  the  shelter 
are  to  be  the  same  as  in  illustration  No.  1  except  that  the  shelter  is 
located  in  Sacramento  instead  of  St.  Louis.  This  means  that  the  out¬ 
side  design  temperature  will  be  94  F  dry  bulb  and  the  wet  bulb  will  be 
69  F. 
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The  sensible  heat  factor  will  remain  unchanged  and  by  following 
the  steps  in  the  previous  solution  the  condition  of  the  entering  air  is 
located  at  point  "H"  (73  db  and  62.1  wb).  Further  calculations  show 
that  the  required  cfm  equals  2550  total  or  5.1  cfm/person,  which  is 
approximately  one-half  that  required  for  St.  Louis.  The  refrigeration 
required  to  cool  the  air  from  the  outside  design  temperature  to  the 
entering  air  condition  at  i,H*'  would  be  4.97  tons. 

In  comparing  the  calculations  for  the  two  cities  note  that  the  out¬ 
side  dry  bulb  temperature  of  94  F  is  common  to  both  cities  so  the  basic 
design  criteria  are  identical  in  the  two  cases  except  for  the  outside 
wet  bulb  temperature.  The  notable  effect  of  this  one  exception 
emphasizes  the  difficulty  in  simplifying  shelter  ventilation  calcula¬ 
tions  which  are  characterized  by  a  multitude  of  variables. 


FUTURE  WORK 

Future  work  on  the  task  will  be  directed  toward  obtaining  more 
data  on  heat  transfer  through  heavy  construction  and  more  information 
on  heat  loss  through  the  floor.  Further  modifications  to  the  psychro- 
metric  chart  are  expected  and  also  some  arguments  on  whether  it  is 
economically  advisable  to  hold  the  inside  temperature  at  85  effective 
temperature  and  90  dry  bulb,  or  at  a  more  comfortable  condition. 
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Anoendlx  A 

*  4 

PERIODIC  HEAT  FLOW  THROUGH  COMPOSITE  WALLS  AND  ROOFS  * 

Alan  Met  Her 


INTRODUCTION 

For  a  given  structure  made  up  of  four  walls  and  a  roof,  each  wall 
and  roof  of  which  is  composed  of  two  or  more  layers,  the  temperature, 
rate  of  heat  transfer  and  time  for  each  inside  surface  can  be  computed. 
In  order  to  calculate  these  variables  various  parameters  must  be  known. 
These  are:  thermal  conductivities,  thicknesses,  specific  heats  and 
densities  of  the  layers,  the  air  film  coefficients  of  heat  transfer, 
the  average  daily  sol-air  temperature  and  the  actual  sol-air  tempera¬ 
ture. 


Three  computer  programs  have  bean  written  which  treat  the  two  layer 
and  three  layer  cases.  One  program  solves  the  problem  for  the  two  layer 
case;  the  remaining  two  programs  solve  the  three  layer  case.  The  three 
layer  case  was  broken  into  two  programs  because  of  the  length  of  the  i 

equations  involved. 

PROGRAM  DESCRIPTION  , 

a.  Program  One  computes  parameters  used  by  Program  Two.  Program 
Two  (three  layer  case)  and  Program  Three  (two  layer  case)  compute  nine 
sets  of  data  for  each  wall  and  roof  on  one  pass,  then  pause,  at  which 
time  the  operator  may  load  a  new  set  of  parameters  and  push  start. 

It  is  desired  to  compute  n  sets  of  data  (n  t  9)  then  one  simply 
changes  the  limits  of  the  Innermost  loop.  Suppose  n  *>  12  then  change 
statement  5+4  from  D071  B  1,9  to  D071  a  1,12.  If  one  would  like  to 
have  data,  say  on  two  walls,  then  the  limits  on  the  outer  loop  should 
be  changed,  l.e.,  statement  4  becomes  D07J  *  1,2. 

Running  time  for  nine  sets  of  data  per  structure  (4  walls  and 
a  roof)  in  the  three  layer  case  is  approximately  1.5  minutes  on  the 
IBM  1620  Computer.  For  the  two  layer  case  the  time  is  approximately 
1  minute.  The  costs  are  about  $1.00  and  $.65,  respectively. 


► 
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b .  SYMBOLS 


Math  Symbols 


t 

e 

Sn 

Si 

t 

0 

A 

e 

<P 

Pi 

p 

m 

po 


Program  Symbols 
CL 
CM 
CO 
HL 
HO 
XKL 
XKM 
XKO 
XLL 
XLM 
XLO 
TE 
TM 
TMM 
TO 

XLAM 

THT 

PHI 

RHOL 

RHOM 

RHOO 
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c.  INPUT 


All  input  variables  are  F  type  format  and  have  a  field  length  of 
10  columns. 

Input  to  Program  1 

Card  1  XKO,  XKL,  XKM,  HO,  HL 

Card  2  CO,  CL,  CM,  RH00,  RHOL,  RHOM 

Card  3  XLO,  XLM,  XLL 

Input  to  Program  2 

Output  from  1  as  is 

Card  5  Tl,  TM 

Card  6  TE,  THT 


Card  14  TE,  THT 
Card  15  Tl,  TM 
Card  16  TE,  THT 


Card  24  TE,  THT 

Input  to  Program  3 

Card  1 
Card  2 
Card  3 
Card  4 
Card  5 


Card  13 
Card  14 
Card  15 


Card  23  TE,  THT 

r  > 


TE,  THT 
Tl,  TM 
TE,  THT 


XKO,  XKL,  XKM,  HO,  HL 

CO,  CL,  CM,  RHOO,  RHOL,  RHOM 

XLO,  XLM,  XLL 

Tl,  TM 

TE,  THT 
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SQUAT! OHS  OF  PROBLEM 

d)EQUAT/OHS  COMMON  TO  BOTH  TWO  LAYER  AMD  THREE  LAYER  CASES 


J  V,  .-  ;  77j  ( 

1  -  J  rt ,  t  l/^-g-^ ; 


o./soe  c. 

"3?  5 

C,  r  CCS  4  COW  7J  y  3/-V7J  S/A//V 


*/ 


c3  s  COS  \  COSH  7rt  -  3/A/  7$  S/NH  7!J  j  Cg  =  3//V7J  COS-V  7J  ;  C+  =  COSTJ  JW  9^ 
Of  s  COS/7  Cafl/V/£V  StN  % S!MH%  /  ^  r  COS  /^C OSH/?  -  3/A/£*  S/MhIq 

a/upcoi*/? ;  <?0 x 

v  s  *° 

*4) ; 

* 8  ft  -  *  *  %  XV*  +c,)-(c,*t> 5  *  jfs,/*  *•  <*)  ; 

4  'ft  -  *4  ijft/j  *  4 1' ft  - » %*»)(<*/%  *4) ; 

Lx)  equATtOMS  FOR  two  layer  case 
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X  = 


^  %  *  *  7c  if}: 

l/a  *  "  l*/aa 

(t.)  «  '  *M  * \  [(*.)  -*«] 

&  r  i'c  6 


t_  =  t, 
m  / 


aeoe(t_  -  tA 

,  y.  \  m _ —J—L 


O.&SS  +  -g- 

«/} 

c)  EQUHTtOA/S  re#  THAEE  LAYS  A  CASE 


W.  Z  COS  Sl„  COSH  fl  f  S/A/  SI  S/A/HSL 
/  *22? 

kV-  =  COS  JV  COSH  Jfl„  -  «W  Jfl  S/A/H  Sl„ 

*  ?  2  ?  ? 

W3  :  S/A/Slf  COSH  Slg  •  s  COS  Slf  3M//V  • 

yH  ^  ^3  ^  ^  s  •%  ^  ^  %  / 


££-  -  7*A/ 


ij^^)*OSf-( 

— * 

(^^)cTcos£- 

("H 

{%  +4 

-  v,SiA  4-  [(^)<&os£-(^^ 

(%*4. 

c/cos£-(^^) 

<6'*£ -**/]}  ’ 
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SAMPLE  INPUT  TO  PROGRAM  1 


1.0 

1.0 

.17 

1.65 

4.0 

.25 

.25 

*193 

140.0 

140.0 

75.0 

.667 

1.67 

.667 

SAMPLE  OUTPUT 

FROM  PROGRAM  1 

.41250 

1.42767 

1*29723 

2.25890 

*1.63005 

2.18189 

*28021 

•  .33777 

22.07677 

*24.51029 

1.42767 

.53511 

2*25890 

*“1.63005 

2.18189 

.28621 

.21986 

20.60052 

2.29368 

11.15752 

5.57533 

.34396 
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SAMPLE  INPUT 


.41256 

1.42767 

.28621 

•  .33777 

2.25890 

41 .63005 

22.29368 

11.15752 

96. 

89. 

82. 

M. 

93. 

-3. 

162. 

-2. 

lid. 

'4. 

114. 

o. 

115* 

i. 

111. 

2. 

164. 

3* 

9$. 

4. 

96* 

93* 

77. 

*4. 

86. 

“3* 

83* 

-2* 

si* 

-1. 

96* 

6. 

116* 

i* 

124. 

2. 

135* 

3. 

i4l* 

4* 

96. 

93. 

126. 

*4. 

125. 

-3. 

lit* 

-2* 

168* 

-i* 

92* 

6. 

93* 

i. 

95* 

2. 

95* 

94* 

4. 

96* 

83* 

7  n 

*4* 

86* 

*3* 

83* 

87* 

TO  PROGRAM  L 
1.29723  2.25896 

22.07677  -?4.Slrt?Q 
2.18189  .28621 

5.57533  .34396 


“1.63065 
1 *42767 
.21986 


2.18189 
.5351 1 
20.60052 
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90- 

0. 

934 

1. 

94. 

2* 

9p« 

f  * 

94* 

I* 

$6* 

160.3 

164* 

*4* 

i  19  * 

*5. 

129* 

13?* 

*1* 

1 42  * 

6* 

144* 

1 46  • 

1  * 

2  * 

132* 

5* 

120* 

4. 

fcp 
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SAMPLE  output  from  program  2 
EXPOSURE  SOUTH 

90.0000uDEGKttS  F .  TM=  fiq./wiMnrcrnrr, 


TE 

THETA 

DGRS  F 

HRS 

82.00000 

-4.00000 

93.00000 

-3.00000 

102.00000 

-2.00000 

110.00000 

-l.oonno 

1 14.00000 

.00000 

115.00000 

1.00000 

111.00000 

2.00000 

104.00000 

3.00000 

99.00000 

4.00000 

EXPOSURE  WEST 

ri=  90. OOOOODEGREES  F. 

TE 

THETA 

DGRS  F 

HRS 

77.00000 

-4.00000 

80.00000 

-3.00000 

83.00000 

-2.00000 

89.00000 

“1.00000 

96.00000 

iOOOOO 

110.00000 

1.00000 

124.00000 

2.00000 

135.00000 

3.00000 

141.00000 

4.00000 

EXPOSURE  EAsf 
TI»  90. OOOOODEGREES  F. 


TE 

Theta 

DGRS  F 

HRS 

126.00000 

-4.00000 

125.00000 

-3.00000 

117.00000 

-2.00000 

108.00000 

-1.00000 

92.00000 

.00000 

93,00000 

1.00000 

95.00000 

2.00000 

95.00000 

3.00000 

94,00000 

4.00000 

TO 

DGRS  F 
89,94869 
89,95005 
89.95116 
89.95215 
89.95265 
89.95277 
89.95228 
89.95141 
89.95079 


Q/A  THETA+PH2/15 
BTU/HR  FT2  HR  MINS 
-.68465  4  32.9 

-.08241  5  32.9 

-.08057  6  32.9 

-.07894  7  32.9 

“.07812  8  32.9 

-.07791  9  32.9 

-.07873  10  32.9 
-.08016  11  32.9 
-.08118  12  32.9 


TM* 


93  »  OOOOODEGREES ' F, 


TO 
DGRS  F 


90.14934 

.24642 

4 

90.14972 

.24703 

5 

90 • 15009 

.24765 

6 

90.15083 

.24887 

7 

90.15170 

.25030 

8 

90. 15343 

.25316 

9 

90.15516 

.25602 

10 

90.15652 

.25827 

11 

90.15727 

.25949 

12 

Q/A  THETA+PH2/1 5 


32.9 

32.9 

32.9 

32.9 

32.9 

32.9 


TMs  93. OOOOODEGREES  F, 


TO 

Q^A 

THETA+PH2/ 

DGRS  F 

BTU/HR  FT2  HR  MINS 

90.15541 

.25643 

4 

32.9 

90.15529 

.25622 

5 

32.9 

90.15430 

.25459 

6 

32.9 

90.15318 

*25215 

7 

32.9 

90.15120 

.24948 

8 

32.9 

90.15132 

.24969 

9 

32.9 

90.15157 

.25010 

10 

32.9 

90.15157 

.25010 

11 

32.9 

90.15145 

.24989 

12 

32.9 

15 
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EXPOSURE  NORTH 

T I =  90.00000DEGREES  F.  TM=  83 . OOOOODEGREES  F. 


TE 

THETA 

TC 

O  y  a 

TU  CT  A  j  pi.''1  ' 

:  .  .  u_  i  i\  i  1  lit-/ 

DGRS  F 

HRb 

DGRS  F 

BTU/HR  FT2 

HR 

MINS 

77.00000 

*4.00000 

89.64615 

-.58384 

4 

32.9 

80.00000 

*3.00000 

89.64652 

-.58323 

5 

32.9 

83,00000 

*2.00000 

89.64689 

*.58261 

6 

32.9 

87.00000 

-1.00000 

89.64739 

-.58180 

7 

32.9 

90.00000 

.00000 

89.64776 

-.58118 

8 

32.9 

93.00000 

1 .OOOUO 

89.64813 

-.58057 

9 

32.9 

94.00000 

2.00000 

89.64825 

*.58037 

10 

32.9 

95.00000 

3.00000 

89.64838 

*.58016 

11 

32.9 

94.00000 

4.00000 

89.64825 

-.58037 

12 

32.9 

HORIZONTAL 

ROOF 

1=  90. OOOOODEGREES 

F.  TM=  100. OOOOODEGREES  F 

• 

TE 

THETA 

TO 

Q/A 

THETA+PH2/ 

DGRS  F 

HRS 

DGRS  F 

BTU/HR  FT2 

HR 

:  MINS 

106  *00000 

-4.00000 

90.50517 

.83353 

4 

32.9 

1 19.00000 

*3.00000 

90.50678 

.83619 

5 

32.9 

129.00000 

*2.00000 

90.50802 

.83823 

6 

32.9 

137.00000 

-1.00000 

90.50901 

.83986 

7 

32.9 

142.00000 

.00000 

90.50963 

.84088 

8 

32.9 

144.00000 

1.00000 

90.50987 

.84129 

9 

32.9 

140.00000 

2.00000 

90.50938 

.84048 

10 

32.? 

132.00000 

3.00000 

90.50839 

.83884 

11 

32.9 

120.00000 

4.00000 

90.50690 

.83639 

12 

32.9 
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SAMPLE  INPUT  TO  PROGRAM  3 

.6 
6.6 
5.33 


1.6 

.17 

.25 

.193 

.667 

6.6 

85. 

64. 

82* 

“4. 

4** 

“*3« 

162* 

si?  4 

1  id  * 

Ul  > 

li4* 

r  * 

li5* 

1* 

1  ii 
164 
44 

65 
77 

66 


83* 

*15* 

84* 

“1* 

46* 

6* 

116* 

1  * 

124* 

2* 

135* 

i* 

l4i* 

4* 

85* 

43* 

I26* 

«4* 

125* 

-3* 

in* 

ii5* 

166* 

*1* 

92* 

6* 

93* 

1  • 

95* 

2* 

95* 

3. 

44* 

4* 

85* 

83* 

1U 

4*4* 

86* 

«3* 

83* 

*2* 

87* 

«*i. 

96* 

0* 

93. 

1* 

1.65 

146.6 


4*0 

75.6 


0.0 
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94* 

2* 

95* 

3. 

94  « 

4. 

65* 

100* 

166  • 

“4* 

119* 

•**3* 

129* 

^2* 

137* 

4* 

142* 

6i 

144* 

1. 

146* 

2* 

132* 

3* 

126* 

4. 
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TI  = 


sample  output 

EXPOSURE  SOUTH 
85 .OOnnnnFaRFEs  F , 


FROM  PROGRAM  3 


TM  = 


TE 

n  i“  r 

THETA 

TO 

UUKO  F 

HRS 

DGRS  F 

82.00000 

“4.00000 

85.15868 

9  3 .00000 

“3.UUUU0 

85.15945 

102.00000 

-2.00000 

85.16007 

110.00000 

-1.00000 

85.16063 

114.00000 

.00000 

85.16091 

115.00000 

1.00000 

85.16098 

111.00000 

2.00000 

85.16070 

104.00000 

3.00000 

85.16021 

99.00000  4.00000 

EXPOSURE  WEST 

85.-15986 

89.00000DEGREES  F, 

Q/A  THET  A+PH2/ 1 5 
BTU/HR  FT2  HR  MINS 
#26182  6  31*9 

.26309  7  31.9 

.26413  8  31.9 

.26505  9  31.9 

.26551  *10  31,9 
.26562  11  31,9 
.26516  12  31.9 
.26436  1  31,9 

.26378  2  31.9 


TI=  85.00000DEGREES  F,  TM»  93 .00000DEGREES ' F, 


TE 
DGRS  F 
77.00000 
80.00000 
83,00000 
89.00000 
96.00000 
110.00000 
124.00000 
135.00000 
141.00000 


T  I  *= 


THETA 
HRS 
-4.00000 
“3.00000 
-2.00000 
“1.00000 
.00000 
1.00000 
2.00000 
3.00000 
4.00000 
EXPOSURE  EAST 
85.00000DEGREES  F. 


TC)  Q/A  THETA+PH2/15 

DGRS  F  BTU/HR  FT2  HR  MINS 
85.31722  .52341  6  31.9 

85.31743  .52376  7  31.9 

85.31764  .52411  8  31.9 

85.31806  .52480  9  31.9 

85.31855  .52560  10  31.9 

85.31952  .52722  11  31.9 

85.32050  .52883  12  31.9 

85,32127  .53010  1  31.9 

85.32169  .53079  2  31.9 


TE 
DGRS  F 
126.00000 
125.00000 
117.00000 
108. 00000 
92.00000 
93.00000 
95.00000 
95.00000 
94.00000 


THETA 

HRS 

-4.00000 

“3.00000 

-2.00000 

“1.00000 

.00000 

1.00000 

2.00000 

3.00000 

4.00000 


TM=  93.00000DEGREES  F, 


TO 
DGRS  F 
85.32064 
85.32057 
85.32001 
85.31939 
85.31827 
85.31834 
85.31848 
85.31848 
85.31841 


Q/A 
BTU/HR  FT2 
.52906 
.52895 
.52803 
. 52699 
.  52515 
.52526 
.52549 
.52549 
.52537 


THETA+PH2/15 
HR  MINS 

6  31.9 

7  31.9 

8  31.9 

9  31.9 

10  31.9 

11  31.9 

12  31.9 

1  31.9 

2  31.9 
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EXPOSURE  NORTH 

TI-  85 • OOOOODEGREES  F.  TM=  83 . OOOOODEGREES  F. 


TE 

THETA 

TO 

Q/A 

THETA+PH2/15 

DGRS  F 

HRS 

DGRS  F 

BTU/HR  FT2 

HR  MINS 

77.00000 

-4.00000 

84.91999 

-.13200 

6  31.9 

80.00000 

-3.00000 

84.92020 

-.13166 

7  31.9 

83.00000 

-2.00000 

84.92041 

-.13131 

8  31.9 

87.00000 

-1.00000 

84.92069 

-.13085 

9  31.9 

90.00000 

.00000 

84.92090 

-.13050 

10  31.9 

93.00000 

1.00000 

84.92111 

-•  130i 6 

11  31.9 

94.00000 

2.00000 

84.92118 

-.13004 

12  31.9 

95.00000 

3.00000 

84.92125 

-.12993 

1  31.9 

94.00000  4.00000 

HORIZONTAL  ROOF 

84.92118 

-.13004 

2  31.9 

I  *  85. OOOOODEGREES  F. 

TM  =  100. OOOOODEGREES' F 

. 

TE 

THETA 

TO 

Q/k 

THETA+PH2/15 

DGRS  F 

HRS 

DGRS  F 

BTU/HR  FT2 

HR  MINS 

106.00000 

-4.00000 

85.59731 

.98556 

6  31.9 

119.00000 

-3.00000 

85.59821 

.98706 

7  31.9 

129.00000 

-2.00000 

85.59891 

.98821 

8  31.9 

137.00000 

-1.00000 

85.59947 

.98913 

9  31.9 

142.00000 

.00000 

85.59982 

.98971 

10  31.9 

144.00000 

1.00000 

85.59996 

.98994 

11  31.9 

140.00000 

2.00000 

85.59968 

.98948 

12  31.9 

132.00000 

3.00000 

85.59912 

.98856 

1  31.9 

120.00000 

4.00000 

85.59828 

.98717 

2  31.9 
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XKG  “  l!i«iuial  couuui.ti'vity  of  cuucieLe,  roomside  layer  BTU/hr  ft  F 

XKL  =  Thermal  conductivity  of  concrete,  outside  layer  BTU/hr  ft  F 

XKM  =>  Thermal  conductivity  of  earth  layer  BTU/hr  ft  F  1 

2 

HO  =  Inside  air  film  coefficient  of  heat  transfer,  BTU/hr  ft  F 

2 

HL  =  Outdoor  air  film  coefficient  of  heat  transfer,  BTU/hr  ft  F 
CO  =  Specific  heat  of  concrete,  rooms ide  layer  BTU/lb  F 
CL  “  Specific  heat  of  concrete,  outside  layer  BTU/lb  F 
CM  =  Specific  heat  of  earth  Layer  BTU/lb  F 

3 

RHOO  **  Density  of  concrete,  roomside  layer  lb/ft 

3 

RHOL  ■  Density  of  concrete,  outside  layer  lb/ft 

3 

RHOM  =>  Density  of  earth  layer  Ib/ft 

XLO  °  Thickness  of  roomside  layer  ft 

XLL  «  Thickness  of  outside  layer  ft 

XLM  *»  Thickness  of  earth  layer  ft 

TI  =  Inside  air  temperature  F 

TM  »  Average  sol-air  temperature 

TE  »  Sol'-air  temperature 

THETA  =  Time,  measured  in  hours  after  noon 

TO  e  Roomside  wall  temperature  F 

2 

Q/A  *  Instantaneous  rate  of  heat  transfer  from  roomside  surface  BTU/hr  ft 

THETA  +  PH  2/15  **  Time  of  day  at  which  TO  and  Q/A  occur  as  a  result  of 
TE  at  time  THETA  (PH2/15  equals  lag  time  in  hours) 
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PROGRAM  1 

READ  101 .XKO  *XKL  *XKM»HO»HL 
READ  iol»CO»Ci-*CM  *RHOO*RHOL»RHOM 
READ  iqi»XtO*xLM*yLL 

H i i “HU /ML 

P! 2*XLO*  it  0 . 1309*RHOO#COyXKO  i**  ( 0 . S  )) 

P 1 3* ( XKO/HO  j  * ( ( 0. ISD^^RHOO^COyXKO) **(0.5) i 

swi«coslpi2i 

SW2»(ExP(Pl2)+EXPi-PI2n*(  .Si 
$W3=sikiPl2S 

SW4*iEXP(PI2  -EXP ( -P I  2 ))* ( .5 ) 

Cl  «iswl*Sw2  +(Swi#$W4} 

£2  ■ ( SWl*SW2  *iSW34$W4) 

<3  »SW3*SW2 
c4  *&Wl4sW4 

GC*U  (£3yPl3)+Cl>**(  2  n  +  <  t  ( C4/PI3 )+C2 ) ** t  2 )  i 

Qc*M2.iAici*Wo.Si 

RC*<C1-M2.*PI3*C4) )*( i C3/PI3 )+Cl i 
RORC+I  (C2-<2.*PI3*C3)  J*UC4^P13)+C2)  i 
SC* I £2*  2.*Pt3*C3i)#i (CS/PlSj+Cli 
SC*SC^I  Cl+<2.*Pl3*C4)i*( |C4^PI3)+C2) ) 
GAM2»XtL*< <O.l309*RHOL*CL>XKLi**i.Sh 
6AM3*<XKL/,HL)*U0.i309*RHOL*CL/XKLl**<  .5)  ) 
Swl»C6s(6AMSi 

SW2«<ExP<GAM2i+EXP(**GAM2))4(,5i 
Sw3*$iM6AM2i  . 

SWJ*>  i  EXP  <  GAMfci  *»EXP  I  ‘■GAM2  i  )# (  .Si 
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D1  »(SW1*SW2)+(SW3*SW4) 

D2  * (5W1*SW2 )—(SW3*SW4) 

03  *SW3*sd2 

K  /.  -  A  i.i  4  li  r  t.J ), 

Mr  .  — '  »i  ^  _'!!  - 

SWl«<£>3/,GAM3)+Dl 

swi=swi*swi 
5W2s(£)4/GAM3  )+D2 
3w2#SW2^$W2 

oo*swusw2 

60*l2./GDi**{.5> 

RD* i D1+ ( 2 .*GAM3*D4 * ft D3/GAM3 ) +D1 ) 
RD-RD+f  <D2-(2.*GAM3#D3>  i*UD4>6AM3  +D2H 
Sb«  ( D2-<  2  .#GAM3#D3  ) )  *  j  l  DS/'GAMS  I  +Dl  I 
St)«SD*-f  <t>l+i2.#GAM3*D4}  i#  UD4>GAM3  Ud2  i  ) 
3  xlks* (XLO/XKO }  +  ( XLLyXKL i  +  j XLM/XKM i 
OMG2*<  (0.i36^RHOM*CM>XKM)*«US)i*XLM 
0M63«IXKM/'H0i*f0.1369*RH0M*CM>XKM|**I.S| 
PUNCH  101 *PI 1 »PI2 »PI 3 »C1 »C2*C3 
PUNCH  16l»C4,GC»RCiSC»GAM2*GAM3 
PUNCH  101 *Dl *D2*D3*d4  *GD*RD 
PUNCH  101»$DiXLkS#OMG2*OM63 
60  TO  1 

101  FOrMAT<6F16.5) 
feND 
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r 


- 


f£y 


1  READ  101*Pil>PI2*Pl3»Cl»C2iC3 
READ  101fC4iQC*RCiSC*GAM2*GAM3 
READ  101 tDl »D2*D3 ,D4tQD>RD 
READ  1 0 1  !OL>')  AU 2  * UrlCj 3 
SWl=(EXP(QMG2)+EXP(-OMG2) )*(,5) 
SW2= ( EXP ( CMG2 ) -EXP ( -0MG2 ))*(.&) 
SW3=S  I N  /  OcMG2  i 
SW4=C0S<0MG2) 

Wl=SW4*SWl+SW3*SW2 

W2»SW4*SW1-SW3*SW2 


W3=SW3#SW1 
W4=SW4*SW2 
PW=(W3*W3 j+(W4*W4| 

PW» ( 2 • /PW ) ## ( • 5 i 

ZW*W2#W3-W1#W4 

YW*W1#W3+W2*W4 

SWla( C3/P 13 )+Cl 

SW2° ( C4/P 13 ) +C2 

ATA*(SWl*SWij+( SW2#SW2 ) 

ATA-(2»/ATA)*#(«5) 

UPS=(C3/PI3)+C1+(C4/PI3 )+C2 

UPS*( (C3/PI3)+C1-(C4/PI3)“C2)/UPS 

UPS*=ATAN  ( UPS ) 

SW4*{PW*PW)*( .5) 

S1a/5»QC#GC 

SW6«0MG3*2. 

SW3*( YW*SW4)+(RC*SW5/SW6) 
SWl®SW3#5W3 

SW7= ( ZW#SW4 ) + ( SC#SW5 / SW6 ) 
SW2«SW7#SW7 

DEL=PW/( (SW1+SW2) ##(  .5) ) 

SWl =SW3 

SW2=(W3+W4)*SW1 

SW4= ( W3-W4 1 *SW7 

SW4«SW2+SW4 

SW5=(W3-W4)*SW1 

SW6= ( W3+W4  J  *SW7 

SW6“£W5“SW6 

SW6-SW6/SW4 
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iBSstJsetMt' 

SW4*5IN(EPS) 

SW3 - P  W k?W 
SW6=QD*QD*( .5 ) 

SWl = ( W3+W4 ) *DEL*SW3 

5W 1  =  SW 1  *-  (  (  W3-W4 )  * DE L* SW4  )  -  Y W 

SW1=P I 1*0MG3*SW5*SW1* ( .5 ) 

SWl=(RD#SW6)-SWl 

swl=swi*swi 

SW2=(W3-W4 j*D£L*SW3 

SW2  =  SW2*-(  (W3+W4)*D£L*SW4i“2W 

SW2*P I l*OMG3*SW5*SW2* ( . 5 ) 

SW2*(SD*SW6 J-SW2 
SW2=SW2#SW2 
SW2=(SW1  +  SW2)*#(  .5) 

GAMA=QD/SW2 
SW7=D3/GAM3+D1 
|  SW8«»D4/GAM3 

SW1=( ( W3-W4 ) *DEL*COS ( EPS! )-( ( W3+W4 )  *DEL*S I N ( EPS ) >-ZW 
SW1=PI1#0MG3*PW*PW*SW1*( .5  I 
SWl  =  (  SD*QD*QD*  (  «  5  )  >-SWl 
SW2=(SW7-SW8-D2i*SWl 

SW3=(W3+W4)*DEl*C0S(tPS)-(W3-W4)*DEL*SIN(EPS>-YW 

SW3  =  P 1 1*0MG3*PW#PW*SW3* ( . 5  j 

SW3  -  (  RD*QD#QD* ( •  5  ) )-SW3 

SW4= ( 5W7+5W8+D2 )  *SW3 

SW5=SW2+SW4 

SW6= ( SW7-SW8-D2 ! *SW3 

SW6=SW6-( (SW7+SW8+D2 ) *SW1 1 

SW6*SW6/SW5 

PSI*ATAN(SW6) 

xlAm=ata*del*gama 

205  PHI=(UPS+EPS+PSI )*(180.)/(3.1416) 

IF (PH  I )6*4*4 
6  PHI=360.+PHI 

4  DO  7  J=1 *5 

5  READ  10 1  *  T I  *  TM 
PUNCH  102»T1»TM 
PUhlCH  106 
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-&■ 


It' 

X. 

ag<- 

;*/ 

V= 

5: 


ni  i  M/~i  i  '  ^  ^ 
r  w  i  »  v.  i  i  X  *-*  -J 

DO  7  1=1*9* 

READ  101 »TE»THT 
SW4=THT+( PHI/15* ) 

if*  l  iit4  ]  vw  I  jx  I 

30  SW4*24 •  +SW4 

31  7MM=T I + ( ( *606* i TM-T I )  )/( ,856+XLKS) ) 

T0-TMM+ ( XLAM* ( TE-TM ) ) 

QDA=1.65*( TO-TI } 

IM=SW4 
SW1  =  I M 
SW2=SW4-SW1 
SW3=(60. >*SW2 
13  IF( IM-12) 10*10*11 
11  IM= IM-12 
GO  TO  13 

10  PUNCH  104*TE*THT  »  T0*QDA* IM*SW3 
7  CONTINUE 
PAUSE 
GO  TO  1 

101  FORMAT  (6F10. 5) 

102  FORMAT (3HTI°F 10 . 5 , 10HDEGREES  F , , 1X3HTM“F 10.5,1 OHDEGREES  F.) 

103  FORMAT (4X6HDGRS  F,7X3HHRS,9X6HDGRS  F,11H  BTU/HR  FT2,3H 
HR,5H  MINS) 

104  FORMAT (2F10 . 5 , 5X2F10 . 5 , 13 ,F5 . 1) 

106  FORMAT (8X2HTE , 5X5HTHETA, 13X2HTO , 7X3HQ /A , 2X1 2HTHETA+PH2 / 1 5 ) 
END 
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PROGRAM  3 

1  READ  101 »XKO.XKL»XKM»HO*HL 

READ  i 01 »C0 »CL  »CM  »  RHOO *RHOL  *  RHOM 

AqUALUj  ALI'DXLL 

PI1=H0/'HL 

PI2=XL0#(  <0.13  09*RHOO*CO/'XKO)**(0.5)  ) 

P  13=  ( XK0/H0  5  *  (  (  0»  1 30 9*RH00 » C0/lXK0  ) *# ( 0 • 5 )  ) 
SWl=C6S(PI2) 

SW2=(EXP(PI2)+EXP<-PI2) )*(.5) 

SW3  =  S I N  (  P 1 2  ) 

SW4=(EXP(PI2)-EXPl-PI2)  )*(.5> 

Cl  =(SWl*SW2}+(SW3*SW4) 

C2  =(SWl*SW2l-(SW3*SW4) 

C3  =SW3#SW2 
C4  =SW1*SW4 

QC~ ( ( ( C3/P 1 3 ) +C1 ) ** ( 2 ) )  +  ( ( (C4/PI3)+C2)#*(2)  ) 
QC= ( <2. )>QC)*#<0.5 S 
RC*|C1+(2.*PI3*C4) )*( (C3/PI3 )+Cl ) 

RC*RC+ ( (C2-<2.*PI3*C3 ) )*l ( CU/Pl 3 ) +C2 ) > 

SC= (C2-( 2. *P I3*C3 ) 1*< (C3/PI31+C1) 

SC*SC~< (C1+(2.*PI3*C4) )*( (C4/PI3I+C2) ) 
GAM2-XLL* ( (0.1309*RHOL*CL/XKL)**(.5) > 
GAM3MXKL/HL)*( ( 0 . 1309*RH0L*CL/XKL ) **< .5) ) 
SW1=C0S ( GAM2 i 

SW2=<EXf?(GAM2)+EXP(“GAM2)  )*( .5) 
sw3*siNiIgaH5) 

SWJ* ( EXP ( GAM 2 ) -EXP { -GAM2 ) !*/ .5) 

D1  =(SWl*SW2)+(SW3*SW4i 
D2  =(SWl*SW2)-(SW3*SW4) 

D3  =5W3*SW2 
D4  =$W1*SW4 
SW1=(D3/GAM3)+D1 

swi*swi*sWi 

SW2=(D4/,GAM3)+D2 
SW2=SW2#SW 2 
QD=SW1+SW2 
QDS i 2 • /QD ) ## ( • 5 ) 

RD»(D1+(2.*GAM3*D4) )*< ( DS/GAMS ! +D1 ) 

RD*RD+  ( (D2-(2.*GAM3*D3) )*( ( D4/GAM3 ) +D2 ) ) 
SD*(D2-<2.*GAM3*D3  )  )*(  (  D3/GAM3  )  +D1 ) 
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SD=SD-( ( Di+(2.*GAM3*D4) )*( < D4/GAM3 1 +D2 1 ) 
2  XLK5={XLG/Xi40)  +  i  XLL/XKL  ) 

SW1= (SD*QD/QC )+( PI 1*SC*QC/GD) 

SW1=SW1*SW1 

SW2=(RD*QD/QC)  +  f  P I 1*RC*QC/QD ) 

SW2=SW?*SW? 

SW2=SW1+SW2 
XLAM=SW2**( .5) 

XLAM= ( 2 • ) /XLAM 

SW1= (SD#OD*OD)+(PI 1*SC*QC*QC) 

SW2=(D4/GAM3)+D2 

SW3=(D3/GAM3)+D1 

SW4= ( RD*QD*QD ) + l P I 1*RC*QC*QC ) 

SW5=( (SW3+SW2)*SW4)+( ( SW3-SW2 ) *SW1 ) 

SW6*(  (SW3-SW2)*SW4i-(  (  SW3+SW2  )  *SW1  1 

SW6*SW6/SW5 

SW6-ATAN  ( SW6  I 

SW1*(C4/PI3)+C2 

SW2* ( C3VP 1 3 ) +C1 

SW3-SW2+SW1 

SW4-SW2-SW1 

SW4=SW4/SW3 

SW4=ATAN ( SW4 I 

203  PHI  * ( SW6+SW4 ) * { 180* 1/(3.1416) 

IF<PHI 16.4.4 
6  PHI-360. +PHI 

4  DO  7  J-1.5 

5  READ  101.TI.TM 
PUNCH  102.TI.TM 
PUNCH  106 
PUNCH  103 

DO  7  1-1.9 
READ  101.TE.THT 
SW4=THT+(PHI/15. 1 
! F ( SW4  j  30.31 .31 

30  SW4*24.+SW4 

31  TMMsT 1+ ( ( . 606* ( TM-T I  1 ) / ( . 856+XLKS )  1 
TO=TMM+(XLAM*( TE-TM) ) 

QDA=1,65*(T0-T I) 

im=sw4 
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SWl= I M 
SW2=SW4-SW1 
SW3=(60. )*SW 2 
13  IF! IM-12J 10,10.11 
11  I M= I M-12 
GO  TO  13 

10  PUNCH  104,TE,THT,TO,QDA.IM,SW3 
7  CONTINUE 
PAUSE 
GO  TO  1 

101  FORMAT (6F10. 5) 

102  FORMAT(3HTI"F10.5, 10HDEGREES  F. ,1X3HTM=F1 0.5,1 OHDEGREES  F.) 

103  FORMAT (4X6 HDGRS  F ,  7X3HHRS , 9X6HDGRS  F,11H  BTU/HR  m,3H 
HR,5H  MINS) 

104  FORMAT(2F10.5,5X2F10.5,13,F5.1) 

106  FORMAT (8X2HTE , 5X5HTHETA, 13X2HTO , 7X3HQ /A , 2X1 2hTHETA+PH2 / 1 5 ) 
END 
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Figure  lb.  Coacret«-«arCh  configurations 


APPENDIX  B  (continued) 
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APPENDIX  £  (continued) 
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APPENDIX  B  (continued) 


CASE  3Zm 


figure  4b.  Concrete-eerth  configuration*. 


Appendix  C 


Table  1C.  Heat  Transfer  Through  Heavy  Wall  Construction 
for  90  f  Inside  Temperature 


Exposure*  8  AM 

10  AM 

12  Noon  ’ 

2  PM 

4  PM 

BTU/hr/sq  ft  « 

Cases  XAt  IIXA,  XVA,  VA, 

VIA,  VXXA, 

VXXXA 

South 

-.0754 

.  A731 

.  niio 

**«•*»• 

weec 

.222 

.223 

.224 

.227 

.229 

East 

.228 

.226 

.224 

.224 

.224 

North 

-.523 

-.522 

-.521 

-.521 

-.521 

Roof 

.747 

.749 

.751 

.750 

.748 

BTU/hr/sq  ft 

-  Cases  IXA, 

XXXB,  XVB,  IVC,  VXXXB 

8outh 

-.155 

-.108 

-.081 

-.088 

-.115 

West 

.378 

.392 

.422 

.487 

.526 

East 

.491 

.471 

.413 

.420 

.417 

North 

-.983 

-.969 

-.953 

-.944 

-.944 

Roof 

1.399 

1.452 

1.486 

1.478 

1.431 

BTU/hr/sq  ft 

-  Cases  XB, 

XXBt  XXXC.  VB 

,  VIB,  VI IB 

,  VXXXC 

South 

-.249 

-.191 

-.158 

-.166 

-.200 

West 

'  .641 

.658 

.695 

.774 

.823 

Bast 

.780 

.754 

.683 

.692 

.689 

North 

-1.618 

•1.601 

-1.581 

•1.570 

-1.570 

Roof 

2.30 

2.36 

2.41 

2.40 

2.34 

BTU/hr/sq  ft 

-  Cases  XC, 

XXC 

8outh 

-.403 

-.158 

-.001 

-.049 

-.196 

West 

.757 

.830 

.989 

1.33 

1.54 

East 

1.35 

1.24 

.940 

.977 

.964 

North 

-2.30 

-2,22 

-2.14 

-2.09 

-2.09 

Roof 

3.25 

3.52 

3.71 

3.66 

3.42 

Worth  Latitude 
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Table  IC.  Heat  Transfer  Through  Heavy  Wall  Construction 
for  90  F  Inside  Temperature  (cont'd) 


Exposure* 

8  AM 

10  AM 

12  Noon 

2  PM 

4  PM 

BTU/hr/sq  ft  - 

Cases  VIC, 

VI1C 

South 

-.342 

-.011 

.187 

.138 

-.061 

West 

.414 

.513 

.728 

1.191 

1.473 

East 

1.225 

1.075 

.662 

.712 

.695 

North 

-1.683 

-1.584 

-1.468 

-1.402 

-1.402 

Roof 

2.36 

2.74 

2.99 

2.92 

2.59 

BTU/hr/sq  ft  - 

Case  VC 

• 

South 

-.617 

.141 

.597 

.483 

.028 

West 

.447 

.675 

1.168 

2.23 

2.88 

East 

2.31 

1.965 

1.016 

1.130 

1.092 

North 

-2.69 

-2.46 

-2.19 

-2.04 

-2.04 

Roof 

3.74 

4.61 

5.18 

5.03 

4.27 

★North  Latitude 


Appendix  C 


Table  IXC.  Heat  Transfer  Through  Heavy  Wall  Construction 
for  85  F  Inside  Teeperature 


Exposure* 

8  AM 

10  AM 

12  Noon 

2  FM 

4  PM 

BZU/hr/sq  ft 

-  Cases  LA, 

XXXA,  XVA,  VA,  VIA,  VXXA, 

VXXXA 

South 

.331 

.336 

.338 

.337 

.335 

West 

.663 

.663 

.666 

.672 

.675 

East 

.673 

.670 

.666 

.666 

.666 

North 

•.168 

-.166 

-.165 

-.164 

-.164 

Hoof 

1.250 

1.254 

1.257 

1.237 

1.253 

BTU/hr/aq  ft 

-  Cases  X1A, 

XXXB,  XVB,  XVC, 

VXXXB 

8outh 

.538 

.584 

.611 

.605 

.377 

West: 

1.071 

1.085 

1.115 

1.180 

1.219 

East 

1.184 

1.164 

1.106 

1.113 

1.110 

North 

-.291 

-.277 

-.261 

-.232 

-.252'' 

Roof 

2.091 

2,143 

2.179 

2.171 

2.124 

BTU/tujtq  ft  • 

Oases  XXXC, 

VIXB 

South 

.371 

.668 

.730 

.711 

.653 

West 

1.133 

1.162 

1.225 

1.360 

1.442 

East 

1.369 

1.330 

1.205 

1.220 

1.213 

North 

-.331 

-.303 

-.269 

-.249 

-.249 

Hoof 

2.298 

2.410 

2.472 

2*481 

$ 

2.366 

BUJ/hr/sq  ft 

-  Cases  XB, 

VS,  VIB,  VXIIC 

South 

.778 

.893 

.967 

.944 

.876 

West 

1.546 

1.380 

1.654 

1.814 

1.912 

Bast 

1.826 

1.774 

1.723 

1.649 

1.643 

North 

_  /.A/. 

■  e  ■*“■* 

-.409 

-.369 

-.342 

-.345 

Roof 

3.104 

3.235 

3.310 

3.321 

3.184 

♦North  Latitude 
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Table  IIC.  Heat  Transfer  Through  Heavy  Wall  Construction 
for  85  F  Inside  Temperature  (cont'd) 


Exposure* 

8  AM 

10  AM 

12  Noon 

2  PM 

4  PM  ' 

BTU/hr/aq  ft  - 

Casas  IC,  IIB,  IIC,  VIIC 

South 

1.184 

1.428 

1.585 

1.538 

1.392 

West 

2.345 

2.417 

2.576 

2.917 

3.125 

East 

2.942 

2.832 

2.527 

2.564 

2.552 

North 

-.708 

-.635 

-.549 

-.498 

-.498 

Roof 

4.835 

5.115 

5.274 

5.298 

5.006 

BTU/hr/sq  ft  - 

Cases  VC,  VIC 

South 

1.139 

1.878 

2.359 

2.211 

1.767 

West 

2.023 

2.425 

2.906 

3.941 

4.569 

East 

4.015 

3.682 

2.758 

2.869 

2.832 

North 

-.921 

-.699 

-.440 

-.290 

-.290 

Rodf 

5.462 

6.313 

6.867 

6.719 

5.980 

★North  Latitude 
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Figure  1.  Graphical  solution  for  equation  (a). 


80 


90 


Dry  bulb  temperature 
Figure  2.  Modified  psychrometric  chart. 


